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C<i ■ Abstract 

. The AjvDY project at RHIC was proposed to measure the analyzing power for 

CN , Drell-Yan production. Test runs took place during polarized proton operations of 

O I RHIC in 2011 and 2012 with a model of the A^vDY apparatus in place. In total, 

an integrated luminosity of 9 pb~^ with beam polarization of 50% was sampled. 



The primary detector components were a hadron calorimeter (HCal) that spanned 
the pseudorapidity interval 2.4 < rj < 4.0 and a small electromagnetic calorimeter 
(ECal). Basic goals for AjyDY test running were to establish the impact of a third 
^ I interaction region on RHIC performance and to demonstrate HCal calibration. 

^ I Energy scale of HCal was established using neutral pion reconstruction and checked 

with hadronic response. In addition, data with a trigger based on HCal energy 
r"! ■ sum were taken to study jet events. First measurements of analyzing power in the 

forward jet production are reported. 

> 

^ '■ 1 Introduction 
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Large transverse single spin asymmetries (SSA) measured in inclusive 
pion production in pp-collisions [H [2] have stimulated significant theory de- 
CN ' velopment to understand the spin structure of the proton. However, this 

process is not easy to describe due to contributions from several mechanisms 
related to initial and final state interactions. Simpler processes from the the- 
^ ! oretical point of view, such as Drell-Yan, prompt photon or jet production, 

I should be considered to disentangle the different mechanisms. In particular, 

inclusive jet production is of interest since it has no final state (Collins effect) 
contribution, and arises only from the Sivers effect. From naive expectations, 
jet SSA should be small because jets integrate over charged and neutral pions, 
and opposite sign asymmetry for vr^ and vr" leads to cancellations. Theoret- 
ical models also expect jet SSA in the forward region to be small [3, 4j. 

First measurements of forward jet production in polarized pp-collisions 
were performed at RHIC 2 o'clock interaction region (IP2) during two A^vDY 
test runs: at ^/i = 500 GeV in 2011 and y/s = 510 GeV in 2012. The A^vDY 
setup in 2011 run is described in [5], and included zero-degree calorimetry 
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(ZDC) for luminosity monitoring and a check of polarization direction. For 
2012 run HCal Left/Right modules were modified into an annular detector of 
20x 12 cells with a 2x2 hole for the beam pipe, and ECal was put on rails so 
that it could be moved out not to shadow HCal for jet measurements. The 
luminosity was 2x higher compared to the 2011 run. Data were taken with 
two basic triggers: 1) energy sum in HCal Left/Right half, excluding two 
outer perimeters to ensure that jets are contained in the detector, with the 
threshold ~ 35 GeV; 2) energy sum in ECal to measure trigger bias on jets. 

2 Jet reconstruction 

The energy scale for HCal was set based on neutral pion calibration [6] 
with a crude adjustment for hadronic compensation from PYTHIA/GEANT 
simulations: E' = 1.12x£' — 0.1 GeV, where E is the incident energy using 7r° 
calibration and E' for HCal cells is used in jet reconstruction. E' = 0.25 GeV 
threshold was applied for a cell to be included into a jet finding. Both ECal 
and HCal cells were used in the jet finding for the 2011 run data. 

Two algorithms were used for jet reconstruction. The cone jet finder starts 
from a seed (high tower in the triggering region of the detector) , sums energy 
in the cone of radius i? = 0.7 in (77— 0) space (where r] is the pseudorapidity 
and is the azimuthal angle relative to the beam direction) around the high 
tower and defines the jet axis from energy- weighted <t]> and <0>. Then an 
iterative procedure is applied until convergence of the jet axis: 1) sum energy 
in the cone of radius R about <riiy>, <(pN> {N is the iteration number); 
and 2) compute energy- weighted <riiy+i>, <07v+i>- 

The recently developed anti-kr algorithm [7j introduces a distance mea- 
sure and uses sequential recombination of cells (clusters) to form jets. For 
cluster pair i and j, the distance is computed as dij = min(A;^j, k^j) x (Rfj/R"^), 
where kT,i = -Ej/cosh(?7j) is the transverse momentum assuming zero mass 
for the incident particle, R^^ = (jji — rjjY + (0j — ^-iid R = 0.7 is used. If 
dij < l/^Tj- for any i, the clusters are merged and the procedure is repeated, 
otherwise cluster j is considered as a jet. Finally, energy and acceptance 
cuts were imposed to select "good" jets: Ejet > 30 GeV, \rijet — 3.25 1 < 0.25, 
\(pjet — 0O//I < 0.5, where (pojf = (vr) for the Left (Right) module. 

Jet shape and transverse momentum distributions from the two algo- 
rithms are shown in Fig{T]for the jet-triggered data and simulations. There is 
good agreement between the data and simulations above the trigger thresh- 
old for both algorithms. There are some quantitative differences between 
the algorithms: cone jets have more steeply falling pj- distribution and more 
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Figure 1: (a) Distribution on jet transverse momentum p-r, normalized by number 
of triggers, (b) Distribution of fraction of energy in the jet as a function of distance 
in {rj—cf)) space from the jet axis (jet shape), (c) Correlation between tower jet 
and particle jet energy. The inset shows the ?7-component of the direction match 
(Ar/) between particle jets and hard-scattered parton, whose direction is defined 
by rjparton, (j^parton- There is 82% match requiring \ A7j\, \ A(j)\ < 0.8. 



narrow jet shape than observed for anti-Zc^ jets. This is hkely related to addi- 
tional "out-of-cone" cells acquired by the anti-fcr algorithm. Hard-scattered 
partons are strongly correlated with jets, and the jet energy scale is checked 
by correlating "tower" jets reconstructed from PYTHIA/GEANT (i.e., full 
detector response) simulations versus jets reconstructed from particles gen- 
erated by PYTHIA, as shown in FigHt. 



3 Jet analyzing power 

The jet asymmetry was calculated from yields in the Left/Right module 
sorted by the polarization direction of the beam heading towards the detector 
for xp > and the opposite beam for x^? < 0: 

e = PBeamA^ = - V^^Vl) / ( + V^^l) , where nI[% 

is the number of jet events in the Left (Right) module for the spin direc- 
tion up (down), as determined from our measured ZDC spin asymmetries. 
This method relies on mirror symmetry in the setup geometry, and cancels 
systematics, such as detector and luminosity asymmetries, through second 
order. To check time-dependent systematics, the asymmetry e was computed 
in xp bins for each RHIC fill and fitted by a constant function, per degree 
of freedom from these fits were close to 1, meaning that systematic errors are 
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small. Bunch shuffling, i.e. random reversing of the spin direction for half 
of the filled bunch crossings that creates effectively unpolarized collisions, 
was used as another estimate of systematics. Mean value of the asymmetry 
calculated for ~100 random patterns was ~(10~^— 10~^), resulting in the 
systematic uncertainty in An less than 2xl0~^. 

The analyzing power was calcu- 
lated using fill-averaged beam polariza- 
tion 0.52 for both beams, and is shown in 
Figj2]for the cone and anti-fcr jet finders. 
The results from the two algorithms are 
consistent within statistical errors. A^ 
at XiT' > is small (~10~^) and positive 
(up to 6a). Jet A^ was also measured for 
ECal-triggered events, and turned out to 
be at the level (1-3)% for > 0. This 
likely means that trigger bias imposed by 
electromagnetic calorimeter prefers jets 
that fragment to a hard neutral pion. 

We have found the jet analyzing 
power to be small and positive. Al- 
though our measurements can help to fur- 
ther constrain the Sivers functions, it is 
most important to measure the analyzing 
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Figure 2: Jet ^a? as a function of 
xp for jet-triggered events. Cone jet 
points are shifted from the mean x p 
value in the bin by -0.01, and anti- 
kp jet points are shifted by 0.01. Er- 
ror bars indicate statistical errors. 



power for Drell-Yan production as a test of present understanding. 
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